u 


u 



!l i 

W ; 

i 

\ I 1 

! I 


1 1 
t i 

y 



|::j 

l -i 

y 



w 


NASA/CR— 1999-208643 


Longwave Radiative Flux Calculations 
in the TOYS Pathfinder Path A Data Set 


Amita Mehta, General Sciences Corp., Beltsville, MD 20705 
Joel Susskind, Goddard Space Flight Center, Greenbelt, MD 20771 



National Aeronautics and 
Space Administration 

Goddard Space Flight Center 
Greenbelt, Maryland 20771 




September 1999 


Available from: 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076-1320 
Price Code: A17 


5285 Port Royal Road 
Springfield, VA 22161 
Price Code: A10 



Abstract 


\ ; 

U 


l i 

y 





H 


!i 

y 




n 



A radiative transfer model developed to calculate outgoing longwave radiation 
(OLR) and downwelling longwave surface flux (DSF) from the Television and Infrared 
Operational Satellite (TIROS) Operational Vertical Sounder (TOVS) Pathfinder Path A 
retrieval products is described. The model covers the spectral range of 2 to 2800 cm- 1 in 
14 medium resolution spectral bands. For each band, transmittances are parameterized as 
a function of temperature, water vapor, and ozone profiles. The form of the band 
transmittance parameterization is a modified version of the approach we use to model 
channel transmittances for the High Resolution Infrared Sounder 2 (HIRS2) instrument. 
We separately derive effective zenith angle for each spectral band such that band- 
averaged radiance calculated at that angle best approximates directionally integrated 
radiance for that band. We develop the transmittance parameterization at these band- 
dependent effective zenith angles to incorporate directional integration of radiances 
required in the calculations of OLR and DSF. The model calculations of OLR and DSF 
are accurate and differ by less than 1% from our line-by-line calculations. Also, the 
model results are within 1% range of other line-by-line calculations provided by the 
Intercomparison of Radiation Codes in Climate Models (ICRCCM) project for clear-sky 
and cloudy conditions. The model is currently used to calculate global, multiyear 
(1985-1998) OLR and DSF from the TOVS Pathfinder Path A Retrievals. 






1.0 Introduction 

‘ ■ !-■, '■ . \ - '• , . ' -4^ • 

Radiative processes are of fundamental importance to die climate system. Outgoing 
longwave radiation (OLR) emitted to space and downwelling longwave flux at the earth's 
Suiftce (DSF) are two of the major components of the overall radiative balance of the earth- 
atmosphere system. It is well recognized that accurate calculations and measurements of these 
components are crucial in understanding their roles in climate variability. 

Several long-term data sets of outgoing longwave radiation are currently available from 
satellite measurements. The longest data set of OLR is based on the Scanning Radiometer (SR) 
and Advanced Very High Resolution Radiometer (AVHRR) measurements taken from 
National Oceanic and Atmospheric Administration (NOAA) polar orbiting satellites. In this 

case, SR or AVHRR measurements of 10 to 12 pm channel radiances have been converted to 

OLR by using regression coefficients pre-determined from radiative transfer calculations 
(Gruber and Krueger, 1984). The NOAA,, OLR data cover the period between 1974 and 
present A similar regression technique to derive OLR from multi-channel High Resolution 
Infrared Sounder (HIRS) measurements has also been developed by Ellipgson et al. (1994). 

An OLR data set has also beat produced from the Earth Radiation Budget (ERB) 
measurements taken aboard Nimbus satellites (Jacobowitz et al. 1984), Hie ERB OLR data 
cover the period between 1975 and 1992 (see Kyle et al., 1993). Another OLR data set has 
been derived from the Earth Radiation Budget Experiment (ERBE) instruments flown cm board 
NOAA-9, NOAA- 10, and the Earth Radiation Budget Satellite (ERBS) (Barkstrora, 1989 and 
references therein). In both the ERB find ERBE data sets broad spectral band radiance 
measurements from die Wide Field of Viqw (WFOV) non-scanner and narrow Reid of View 
(NFOV) scanner instruments have been used to derive two sets of OLR. The ERBE science 
team also provides estimates of clear-sky y^ucs qf outgoing longwave radiation (COLR) based 
on tile scanner measurements (Barfcstropi etaL 1989). The ERBE scanner OLR and COLR 
data are considered to be of very gopd quality in terms of their spatial and temporal coverage 
and instrument accuracy. But die global ERBE scanner data are available only for die limited 
time period from February 1985 toMay LiSd. , 

OLR, as menticmed above, has been obtained from top-of-atmosphere satellite 
measurements. DSF, on the other hand, has no significant correlation with the top-of- 
atmosphere radiation (Ramanathan, 19$6), ,and so far has not been derived directly from 
satellite measurements. Also, surface-based measurements of global DSF are not possible. A 
plausible approach to obtain global, logg^aim DSF is to calculate the DSF from satellite- 
derived temperature, moisture and cloud information in radiative transfer. For example, Gupta 
et al. (1992) have developed a parameterization to calculate DSF based on radiative transfer 


from tovs, 

and clouds from tbe ly ^n ^ on al Cloud Qyr^o gy Pipject (ISCCP). Darnell et al. (1992) 
have used this parameterization to calculate o$^y DSF fa 

The present study is an attempt to produce global, long-tom records of OLR and DSF 
based on satellite measurements. As described by Susskind et al. (1997), global multi-year 


(1985 to 1998) PathfinderPath A retrieval products doived from the TOVS radiances are 
currently available. We , developed a radiative transfer model to calculate OLR and DSF as a 
function of relevant geophysical parameters in producing the longwave radiative fluxes in the 
Pathfinder Path A data set Using the model and the Pith A retrievals, global, twice-daily per 
satellite fields of OLR, clear sky OLR (COLR), and DSF are calculated and included as a part 
of the Pathfin der Path A data set This approach is similar to that of Wu and Susskind (1990), 
but a new, computationally efficient radiative transfer model with improved accuracy has been 
developed. Wu andSUsSkind used monthly mean soimdings to calculate monthly mean OLR 
based on TOVS retdmriis, wlnfe the Path A OLR, clear dry OLR, and DSF have been 
on a sounding by sounding basil Also, Wu and Susskind used a climatological 
ozone distribution, whereaswe use die TOVS ozone retrievals in the calculations. 

The objective of the present rcport is tb describe the radiative transfer model that we 
have developed for the calculation of OLR and DSF using the Path A retrievals. It is well 
known that catenations of OLR and DSF require integration of radiances arising from 
numerous spectral hoes; and from all zemth and azimuthal directions. Therefore, radiative 
transfer models are commonly used to estimiteOLRaiid DSF. Most radiative transfer models 
are developed to calculate longwave flux for spectral bands rather than for monochromatic, 
individual spectral lines. This is accomplished by modelhig the atmospheric 

transmittances for spectral bands. Examples of these are statistical band models (e.g. Briegleb, 
1992 and references therein), the correlation k-distribution method (e.g. Lacis and Oinas, 


1991), and pre-computed t ransmittanc e table lbbk-up methods (see Chou aid Kouvaris, 1991 
and references therein). The radiative transfer model we developed is a medium spectral 
resolution band model. In our model, the param eterization used to compute band averaged 
t ransmittan ces as a function of temperature, moisture, and ozone profiles is a refinement of the 


approach originally used by Susskind et al. (1983) to model channel transmittances for the 
HIRS2 instrument of TTROS-N. Also, the parameterization methodology is analogous to that 
used for nwtftHng the TOVS channel transmittance in the analysis of TOVS data in the 
Pathfinder Path A data set 

In addition to the frequency integration, calculations of OLR and DSF require 
(hrectional integration. Directional integration of transmittances or spectral radiances is 
generally carried out by the diffusivity approximation. According to this approximation, it is 
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assumed that angle-averaged, band-averaged transmittances can be accurately represented by 
the band-averaged transmittance of a beam along an effective zenith angle 9 =52.96°, where 
sec 0=1.66 is the well known diffusivity factor. The diffusivity factor was proposed by 
Elsasser (1942) and has since been used in radiative transfer calculations for meteorological 
applications (see Goody and Yung, 1989 for a discussion on this topic). In band models 
where band transmittances are generally parameterized in terms of optical mass (effective 
absorber amount), this is accomplished by simply multiplying the band optical path by the 
factor 1.66 in the radiative transfer calculations (Stephens, 1984). In our model we define and 
obtain an effective zenith angle for each band used in the flux calculation such that the band 
radiance in that direction, multiplied by 7t, is essentially the same as the radiance integrated 
over all angles. The use of effective zenith angles can be thought of as a band-dependent 
diffusivity approximation. 

In this report, we describe the radiative transfer model, including the methodology used 
to obtain the effective zenith angles. Accuracy of the model used to calculate OLR and DSF is 
assessed by a comparison with line by line calculations. Moreover, our model calculations of 
OLR and DSF are compared with a number of other radiative transfer model calculations 
provided by the Inter-Comparison of Radiative Codes for Climate Modeling (ICRCCM) 
program (Ellingson et al., 1991). Comparisons of the TOVS Path A OLR with ERBE and 
AVHRR OLR, and that of Clear-sky OLR with ERBE clear-sky OLR have been presented by 
Mehta and Susskind (1999). Comparison of the Path A DSF with available surface 
measurements is underway and will be reported elsewhere. 

The radiative transfer model used to calculate OLR and DSF is outlined in section 2. 
Comparisons of model calculations for various cases with line by line OLR and DSF 
calculations and with the ICRCCM results are provided in section 3. A summary is given in 
section 4. 


2.0 Overview of the Calculation of Outgoing Longwave And 
Downwelling Surface Radiation 

OLR is the upward thermally emitted flux, pt , at the top-of-atmospherc and DSF is the 
downward thermally emitted flux, F^ , at the surface. These are related to monochromatic 
radiances according to: 

^ 2tc jc/2 «> 

OLR= F ' = J J J R v ~ 0,4>) sin0 oos6d6tydv 2.1a 

4> = 0 0 = 0 v»0 


3 



and 


I 2 n a/2 °o | 

DSF = F 4 '= J J J RJ (0, 

f»09=0v=0 


2.1b 


where 6 and <(> r epr esent zenith and azimuthal angles respectively; R; is the thermally 
matted m »wchBw ^fe radiance going to space, and R* Is downwelling radiance at the 

surface at these angles and at wavenumber v . Rot a cloud-free atmosphere, with unit surface 

t a* 

emissivity, outgoing radiance R v can be expressed as: 


2.2a 


hT* (9 , ♦) = B w (T f ) < (P. . 9 , ♦) + )b v [T(P)] dT »(^’ t ? dP 

p * 

Fwactoudy,<wen»st«diiospherew&speairidoodettrf^ Rl 041 can be expressed 


Rf*(t„A*)=e. 


B,(T t )t:(P I ,8,4,)+) B.[T(P)1 jpj + (t-.a.)R^“(9 . +) 

2.2b 


p. 


if it is assumbd that the clouds have a trananisrfvify given by (1-ecv )• Similarly, downward 
radiances for a cloud-free atmosphere at die surface, Rj 0 * , can be expressed as: 


Rf* (9 , ♦) = I B.[T(P)] dP 


2.2c 


and fear overcast atmosphere as: 

(e CT .e 


, ♦)=e„[B.(T,K(P.. ».♦) + ! B,[T(P)] : ^^fldpl + (l-e.,)R;"(9,+) 

. p* 


2.2d 


In die above equations, T g is the surface temperature, T c fte cloud-top temperature, T b 
the cloud-base temperature. By die Planck function evaluated at wavantmber v, P t an effective 
pressure of the top of etmo*$>here, Pc and P b the cloud-top and cloud-base pressure, 
xj(p ,6 ,$) is die atmospheric transmittance from pressure P to the top of the atmosphere. 
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and Xy(P,6,4») is the transmittance from pressure P to the surface. (Also, note that 
everywhere in the text, upward pointing arrows represent outgoing radiance or flux at the top- 
of-atmosphere, and downward pointing arrows represent downward radiance or flux at the 
sfirface). 

If for a given scene, the cloudy portion of the scene has fractional coverage a, and die 
clear portion 1- a, then upward and downward radiances for the scene can be expressed as: 

R ( v TorA) = (1 - oe^jRr” ;) (0 , + as^ ”^(0 , if* 2.3 

where R(^° r ^^(0,<|)) CLD are the upwelling and downwelling radiances which would arise 
from the scene if it were covered by a cloud with emissivity of 1, and die parameter ae c is the 

radiatively effective cloud fraction. Combining Equations 2. 1 , 2.2, and 2.3, and ignoring the 
frequency dependence of , gives: 

OLR = (l-ae c )F TCLR + ae c F TCLD 2.4a 

DSF = (1 - at c )F l CLR + ae c F A CLD 2.4b 

>.CLR | CLR 

where F 1 and f 4 ' are the top-of-atmosphere upward flux and downward surface flux 

a Pt.ti I pt.ti 

from the cloud-free portion of the scene, and F and F^ are the analogous fluxes 

which would emerge from the cloudy portion of the scene if the cloud emissivity were 1. 

According to the above set of equations, OLR and DSF calculations require integration 
of radiances over angles 0, and <j> , and over wavenumber v . The radiative transfer model we 

developed accounts for the angle and frequency integrations in the OLR and DSF calculations 
using the following approximations : 

• The radiances are considered to be homogeneous in the azimuthal direction <)> . They 

are dependent only on zenith angle 0. 1 f 

• Integration over the zenith angle is approximated by evaluating radiances within a 
spectral band at single band dependent effective zenith angle, and multiplying by jc to 
obtain the flux. The effective zenith angle is found from an ensemble of atmospheric 
profiles such that tire angle integrated flux best approximates what would have been 
obtained if radiances were integrated over all angles (see section 2.2). The effective 
band zenith angles are found separately for OLR and DSF calculations. 

■ . : - . - , - - , r ' "f . ‘ 
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• The spectral range considered for the calculations oxfends from 2 to 2750 em~l . The 
frequency integration is carried out as a summation of fluxes over 14 spectral intervals 
or hands, with band centers at frequency v, Wm, each interval, representative effective 
zenith angles 6, for OLR, and 0* for DSF calculations are found based on 
monochromatic radiance calculations. Upward and downward transmittances for each 

at the appropriate effective zenith angle for that interval using a rapid algorithm band 
transmittance model. In die rapid algorithm, band transmittances and ^ are 
parameterized as a function of temperature profile, water vapor profile, and ozone 
profile (see Section ..23). 

• The Planck function for a spectral interval is approximated by its value at the center 
wavenumber of that band, Vj , such that Bi(T)-Bi(Vi,T). 

Using these approximations, clear-sky OLR and DSF are evaluated as : 


ptCLR _ T7 F tajt = R £ Bi(T s )tf(p 8 ,0f)+fBi| 
i=l i=l », 


,dxt p.Of 
1 < dP 


F* CLR =iFj CLR = ic £j Bi(TXP))^^dP 


where tj|p , 0^) are the atmospheric transmittances evaluated W band-averaged effective 

t _-:£i . ■=: f«,;! -yoi> 

zenith angles 0, and 0, . Analogous equations can be written for F 1 and F* . An 

t , CLD 

assump tion of unit cloud emissivity in the calculation of F 1 W is made because the cloud 

emissivity is included in the frw <XEc in equations 2.4a-b. 

The effective zenith angiesand band transmittances are derived from monochromatic, 
line by line transmittance calculations. In tie following subsections, the line by line 
rnimitt rinn g, method of finding the effective zenith angles, and the parameterization of the bold 
transmittances are briefly described. - 


2.1 Line by tine Calculations 

Monochromatic line by lie (LBL) 
carried out for a reference climatological i 


sorption coefficients used in this study were 
If' of temperature, water vapor, and ozone, 
at 40°N. the muring ratios of CO 2 , N 2 Q, and 
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CH4 were assumed to be constant with height in this profile, with values given in McClatchey 
et al. (1972). For the reference profile, LBL calculations of integrated optical depths in 66 
atmospheric layers from 0 to 1000 mb were obtained between 0 to 2750 cm' 1 at every 0.01 
cm' 1 (Susskind and Searl, 1978). Optical depths for two additional temperature profiles, one 
colder than the reference temperature profile by 20°C and another warmer by die same amount, 
were also computed. LBL optical depths for any temperature profile were modeled as a 
quadratic function of temperature based on values calculated for the reference profile and the 
two additional profiles. Monochromatic optical depths for any water vapor (ozone) profile 
were computed by multiplying the reference optical depths values by the ratio of the layer water 
vapor (ozone) column density in the profile to that in die reference profile. The total 
transmittance function also includes absorption by the Nitrogen continuum and water vapor 
continuum, as in Susskind and Searl (1978). Monochromatic optical depths at any zenith angle 
were obtained from the nadir values by multiplying them by secant of the zenith angle. 

In order to generate the effective zenith angles and the coefficients of the band 
transmittance model, 22 different profiles of varying temperature, water vapor, and ozone were 
used. Monochromatic transmittances (and radiances) for these profiles at eight zenith angles 
were calculated as described above. The eight angles chosen for die calculations, have been 
used in the angular integration of radiances by Gaussian quadrature. 


2.2 Effective Zenith Angles 

To obtain the effective zenith angles for model bands, the concept of monochromatic 
effective zenith angle is first defined. As mentioned earlier, azimuthal symmetry of radiances is 
assumed. For a given atmospheric profile n, if the outgoing and downwelling surface 
radiances were independent of zenith angle and were equivalent to their values at angles 
6e <V' n ) and de respectively then: ; : ’ ,- :i 



The single profile effective zenith angles, tijn) and 8e (v#n) are therefore defined so that 
the radiances and calculated at these angles and die actual angle-integrated radiances 
F v t and Fy satisfy Equations 2.6a-b. Radiances generally vary monotonically with zenith 



angle for upwelling (downwelling) flux if the 
monochromatic radiance b coming from the troposphere. In fee e ase of upwelling radiance, 
fte reverse is true if the radiance is coming from the stratosphere. Consequently, the effective 
zenith angle ft my wavenumber hi well defined for a given profile. As described in section 
2.1, line by line monochromatic upwelling and downwelling radiances Rj and for eight 
Gaussian quadrature angles and corresponding fluxes, Fj and F* have been computed at a 

0.01 cm -1 spectral spacing for 22 atmospheric profiles; Using these calculations, effective 
zenith angles, 0*( A )(v,n) > ^ or each wavenumber and profile were found by linear interpolation 


of R(6) in cosme(O), so that Rj equals and rJ equals for each case. The 
values of effective zenith 0[ and Of for the model spectral band i arc computed according to: 


Vi 22 Qft’i 

£ £eJW(v,n) — * 


dR T W(v,n) 


d6 


fl T(i) = v = v. n = 1 

e ’ V* 22 dR^Cv.n) 


2.7 


Vi 

t X 


de 


v = v. n = 1 

dRtW(V4t) 

where vl and v 2 are the lower and upper limits of wave numbers for bandi and 1 --V — 

CB 


is the derivative of profile n radiance with respect to angte at ^ dfrctiye zphh angle at each 
frequency. Equation 2.7 ensures that the profiles and frequencies for which radiances vary 

most wife zenith angles cany appropriately more weight. 

Figure la shows monochromatic spectral outgoing flux Fj , and Figure lb shows the 

equivalent monochromatic brightness temperature spectrum, given by 0 v »Bv(”“)» 


averaged over the values obtained from the 22 profiles, in the spectral range of 2 to 2750 cm -1 
Analogous spectra for the downwelling surface flux F^ , and corresponding equivalent 


monochromatic brightness t em p e ngnie are shown in Figures lc and Id. The upweiliQg flux 
spectra given in figures la and lb deafly show absorption bands of water vapor, carbon 
dioxide and ozone, along with atmospheric windows in regions between 800 and 1000 cm'* 
and between 1050 cm' 1 and 1200 cm' 1 . The downwelling flux spectra in Figure lc and Id 
show features which are out of phase with those of the upwelling flux spectrum. As a result of 


the strong absorption by OO 2 and water vapor near the earth's surface the emission level is 
pl oa * to the surfaoe over most of the spectrum, except in the atmospheric window regions in 
which the atmosphere is not opaque and some radiative fluxes coming from colder 
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Figure 2 shows monochromatic radiances rJ (9) and rJ (9), averaged over die 22 
profiles evaluated at each of five different zenith angles in the same spectral range. These 
- angles are a subset of angles used in die angular integration of radiances by Gaussian 
’ quadrature. The outgoing radiances (Figure 2 top panel) generally decrease with increasing 
zenith angle. This is due to increased attenuation of radiation arising from warmer, lower 
layers in the atmosphere and increased emission from colder, higherlayers asr angle increases, 
f; However, at the wavenumbers where the main emission going to space arises from die 
stratosphere, e.g. in the center portion of 15 pm CO 2 band, the radiances increase with 
! ' increasing zenith angle. The downwelling radiances, on the other hand, ate independent of 
!; tztatftlt angle over most of the spectrum but increase with increasing zenith angle in the 
atmospheric window region (bottom panel of Figure 2), because as the optics! depth increases 
Mgfcfef zenith angles, more downwelling radiation arising from the lower, wanner 

To give an indication of the corresponding monochromatic effective zenith angles, we 
define and plot (v) , such that: 


22 


• C X®! W (v,n) 


dR T W(v,n) 


dO 

g dR;W(v,n) 
n = 1 


de 



2.8 



These aisles are shown in Figure 3. The solid horizontal Ime in Figure 3 is # 52.96°, which is 
the frequency independent angle commonly used in the diffusivity approximation. At the 
frequencies for which die derivative of the radiance with angle is less than 
l.Oe* 5 mW»m‘ 2 »ster‘ 1 *(cm' 1 )* 1 •degree' 1 , die effective angles are not shown in the figure since 
radiances arc almost the same at all angles in these cases and do not contribute to the band- 


averaged effective zenith angles; As shown in Figure 3, effective zenith ingles for outgoing 
flux ® e (Y) , in general, remain between 48° to 58° in die large portion of the range 2 to 2750 
cm'l. In the window region between 800 to 1250 cm~l, the angles remain approximately 
between 55° to 58°, except in the ozone absorption band where they are between 40° to 52°. 
Effective angles in regions of strong water vapor absorption tend to be about 50°, and values in 
die 15 pm and 4.3 pm CO 2 bands are somewhat lower. The effective zenith angles for 


downwelling flux, ®g(v), are lower compared to those for outgoing flux. In the case of 
downwelling flux, emission is primarily coming from the bottom of die atmosphere where lines 
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are stronger. In 


regimes in which radiances vary considerably with zenith angles, 
e.g. between 8QQ and 1200 (mi- 1 , values of ?eM are found to be between 52° to 56°. In the 
remaining part of the spectrum, where the variasien of downwelling radiances with angle is 
negligible (see Hgure 2), effective angles are between ^ to 45°. The ttuyor points of interest 
are dun the effective zenith angle varies considerably with frequency, the values differ from 


The behavior of effective angles can be imderstood based on how radiances vary with 
Inequation 2.1, if radiances R v (0) were linear is seed, it is easy to show that 
the 6 e is equal to 60°. This would occur if toe monochromatic absorption coefficient was 

small raid toe atmospheric transmittance can be approximated as linear in the optical path. If 
R v (0) were linear in 6 (not a physical situation), then 6 e would be 45°. To the extent that 
r v (0) varies faster feat linear in 0, 0g would be between 45° and 60°. On die other hand, 
0 e < 45° would occur if R v (0) varied slower than linear in 6. This can occur for very large 
absorption features which saturate at intermediate zenith angles. 

The 14 bands chosen for the calculations, and the means and standard deviations of toe 
effective zenith angles for these bands based on the 22 profiles, are given in Table la for OLR 
and in lb for DSF, along with corresponding values of secant(0i). As mentioned earlier, for 
computational purposes, the frequency-independent value of secant(0)=1.66 is generally used 
in most radiative transfer models. We find that toe bend values of secant(0;) for OLR range 
from 1.57 to 1.76, and the average value is close to 1.66. However, for DSF many values of 
secant(0i) are significantly lower than 1.66. 

Using toe effective zenith angles for each braid shown in Table 1, band-averaged 
UArahittinces are computed for toese angles asa fraction of the teragienstrae, moisture, and 

section. 1 : ■ r "-' ; v ' 


2.3 Algorithm Band Transmittance Model 



used by Susskind et al (1983) to 


ftROS-N. Analogous methodology for modeling the TOYS channel transmittance is used in 
analyzing the TOVS data in toe Pathfinder Path A data set (Susskind et al 1997). 

As in Susgkmd et al (1983), we select a form for band i averaged transmittanoes t,(P), 
which reflects the behavior of monochromatic transmittances x v (P), where t,(P) represents 
the t ransmittance from pressure P to the boundary pressure P b , and P b is the top <rf toe 
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1 


atmosphere in the case of upwelling radiation, and the surface in the case of downwelling 
radiation. The monochromatic expressions are exact and given by the two basic relations : 

*.(pj)=rw p i- p H)' 29 

j-i 


where t v (Pj , Pj_i) is the transmittance in a layer between Pj and Pj-.* and 


T .( P J ■ p h) = <( p J • p hK ( p i- p h) <( P i • p h) C( p i . p n) <?(*, • P J-.) 210 


where x f , x w , and x° represent transmittance of gases with fixed distribution, water vapor, 

and ozone respectively, and x** and x" 6 represent transmittances resulting from water vapor 
continuum and nitrogen continuum absorption. 

Our goal is to model the band averaged transmittance x,(P) to be used in Equation 2.5 
for the computation of fluxes as a function of atmospheric temperature and constituent profile. 
It is convenient to use equations of the form of equations 2.9 and 2.10. But equations of the 
form of 2.9 and 2.10 do not hold for band averaged transmittances unless the transmittances 
under question are constant across the band, or the terms being multiplied are uncorrelated in 
frequency. We therefore define effective band i transmittances, 

Xi(Pj , Pj_,) , x^(Pj, P H ) x* (Pj , Pj_, ) , and x?(Pj , P H ), in a manner such that the analogous 
equations hold: 


*,( p ) = rii'( p j. p j-.) 211 


i.(P, . P H ) = (P, . P H ) . ^ (P,.P H ) if (P, . P H ) . 2.12 

Note that continuum absoiption varies only slowly with fiequeney in a band so the inclusion of 
band average values for those terms in eqttatipns of the form of 2.9 and 2.10 will not introduce 
an appreciable error in the use of analogous expressions for band averaged transmittances. 

To ensure the form of Equations 2.1 1 and 2.12, we define 
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^(P, , P,-.) - <**1**,K 2.13 



2.15 




in 


theiayer Pj , P H . For band averaged effective transmittances, all forms; of ii(p J t P^,) also 
depend on conditions between P H and the boundary pressure P b . The coefficients of the model 
are generated based on a set of line by line calculations, in which all terms in Equations 2.13-2.16 
are computed for the ensemble of 22 profiles, each done at the appropriate zenith angle 0, for band 
i (as described in Section 2.2). 

The form of the rapid algorithm parameterizes xl (Pj , Pj_, ) , x* (Pj . Pj_i ) . and x?(Pj ,Pj_,) 
in terms of the mean temperature in layer Pj , Pj-p the column density of water vapor and the 
column density of 0 3 in the layer P } , P^,, and accounts for die temperature-water vapor-ozone 
profile between P J _,and P b . Appropriate coefficients are generated for both OLR transmittances. 


for which P k = 0 mb, and DSF transmittances, for which P b = 1050 mb. Note also that the 
defitddonof t^Pj.Pj.,) ^different for OLR and DSF. For all terms, weparameterize the 
effective layer optical depth, In Xi ( p j * y- as a linear combination of the variables. 


12 


\ ! 

U 


! 

W 



U 

y 


0 

\ h \ 

y 




u 

M 






In the case of fixed gas absorption, the transmittances depend only on the temperature 




profile. We model 

to if (Pj , Pj-,) = Ay+ By (ATj, H ) + Cy (ATy.,) 2 + DyAT, 


ABOVE j-1 


2.17 


ATj ^j in Equation 2.17 is given by the difference between (be profile temperature in the layer 

between Pj and Pj rI , and that of a standard temperature profile T®j_,. AT above ._ 1 in 

Equation 2. 17 is a weighted difference of the temperature profile above (below for downward 
flux) Pj_, from that of a standard temperature profile 


AT 


ABOVE j-l 


H AT k , k - 1 (t[(P k )-xf(P k _ 1 )) 
k=l 1 — X| (P k _i) 


2.18 


which weights layers above (below for DSF) Pi- according to their importance in die radiance 
integral. For each band, and layer, we have 22 linear equations to determine the four 
coefficients. The coefficient D y is set equal to zero forj=l, and for cases in which 3 ft)* 

greater than 0.99, indicating that the levels above j do not contribute significantly to die 
radiance. • •;? 


For the water vapor and ozone, we parameterize the effective absorption coefficients in 
terms of the water vapor profile and ozone profile respectively. The model at this time ignores 
the temperature dependence of die water vapor and ozone line absorption. As in the case of die 
the fixed gas transmittance (equation 2.17), the profile is referenced to the standard water vapor 


and ozone profiles u^ ^, and uV_, according to ij 

layer density in the given profile. 

For water vapor we write - 



where u Jtj _, 


is the column 


Inti (Pj »Pj-i) — Ey r w j + I^j r w j + Gy r w j + Hy r w above , j-l 2.19 


where r wABO VE,j-i is a weighted esthwite of the effect of die water vapor ratio above (below 
for DSF) Pj_j, given by 
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r wABOVE,j-l 


2 ® 


2.20 


where iL (P H ) is the effective water vapor transmittance for the standard profile. As with 
temperature, we have 22 equations to determine four unknowns in each layer. Hy is set equal 
tozeroifj- 1 or tLfPj.,)^ 0M The ozone algorfibih is identical in form to the water vapor 
algo rithm, with coefficients I, J,K, and L for band land layer j,j-l, except that coefficient L 
multiplies a slightly different term for r oABOVEj _, given by 


r o ABOVE, j-1 — 



2.21 


The form of this term for ozone was found to fit the line by line cal c ulations better than die 
form used for water. 


Based on die parameterized transmittances, dear sky OUR and DSF are obtained 
according to Equation 2.5, with analogous equations for dbudy OLR and DSF. In the 
following section,the accuracy of die model calculations of OLR and DSF are presented by 
comparing diem with die line by line calculations carried out by us and also with those 
provided by the 1CRCCM results. 


3.0 Accuracy of the koagwave Fluxes 

3.1 Comparison with line by line Ca lc ula tio ns 

Given the accuracy of the line by line calculations, the calculations of OLR and DSF are 
based on two approximations: use of die effective zenith angle for each band to account for 
angular integration; and use of the rapid algorithm for band transmittance, together with the 
evaluation of the Planck function for each band at its central frequency, to account for 
frequency integration. If both approximations were perfect, band fluxes would agree exactly 
with those obtained from lime by line calculations fra- the same profiles* It is therefore 
necessary to assess how well die model calculations compare with die line by line calculations. 
This c om p ari son should be done for both clear and cloudy conditions. For the cloudy 
comparison, a cloud with emissivity=l is placed at the 500 mb level in each of the 22 
profiles. OLR and DSF are then calculated for clear and cloudy conditions by the tine by line 
method and by using the model described in the previous section. 
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Results of the line by line and the model calculations 6f OLR and DSF for the clear-sky 
and cloudy cases for the 22 profiles are shown in Figure 4. Corresponding means and 
standard deviations are given in Tables 2a and 2b for both total OLR and its spectral 
components, and in Tables 2c and 2d for DSF. As shown in Figure 4, out model calculations 
of OLR(left panel) and DSF (right panel) for clear-sky and cloudy conditions agree well with 
toe corresponding line by line calculations. As shown in the Tables 2a, 2b, 2c, mid 2d, 
agreemeritfor all toe spectral components is also good. The mean differences between toe line 
by line and the model calculations of total OLR under dear conditions is 0.26 W*m-2 and the 
standard deviation is 1.31 W»m-2. For the cloudy conditions, toe mean difference of total 
OLR is 044 Wnn-2 and the standard deviation is 1.89 Wnn-2. For total DSF, toe mean 
differences arc -0.71 W*m-2 and -0.3 W*m-2 for clear and cloudy conditions respectively. 
Corresponding standard deviations of DSF arc 0.56 W*m-2 and 0.18 W»m*2. These errors 
result from a combination of errors due to errors in the frequency and angle integral 
approximations, as well as the band transmittance algorithm error. ^ : ~ ’ < 

As described in the previous section, we used toe band-dependent effective zenith 
angles to approximate toe angular integration required to calculate OLR and DSF. The 
conventional diffusivity approximation uses a frequency-independent effective zenith angle of 
52.96°. Band coefficients are found which represent nadir transmittances (effective zenith 
angle of 0°) and when used in toe calculations the nadir optical path is scaled by, the factor 1.66 
=s secant (52.96°). As described in Section 2.3, the model we use ; to calculate toe band 
transmittance was originally developed and used for H1RS channel transmittances. In the 
model, layer band transmittances for water vapor and ozone are parameterized as a non-linear 
function of toe ratio of layer column density to toe column density of toe tgiprefiriate gas in the 
reference profile. Therefore, in the context of our model, toe conventional diffusivity 
approximation, in terms of scaling the optical path, is not directly applicable. Instead, line by 
line calculations for each band arc carried out at toe appropriate effective zenith angles to obtain 
the band transmittances. 

To examine different approaches of angular integration, we carried out two additional 
sets of OLR and DSF calculations. In toe first case, the band dependence of toe effective 
zenith angle is ignored, ami toe monochromatic transmittances arc calculated at a constant 
zenith angle of 52.96° and then averaged for use in toe parameterization of the band 
transmittances. The second case uses a different approach, in which monochromatic 
transmittances for a given profile are first integrated over all zenith angles using Gaussian 
quadrature and tom averaged to obtain toe band transmittances. In both these cases, 
appropriate rapid algorithm model coefficients based on toe line by line calculations, as 
described in Section 2.3, are generated. Total flux values computed by using each of these 
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approaches, along with those obtained from line by line calcu l a ti o ns aid our own model 
calculation* (as shown in Table 2a and 2c) are given in Table 3. As shown in Table 3, in the 
first additional case (column C), the mean difference of OLR from line by line calculations is 
-0.61 Wnnr 2 (standard deviation of 1.63 Wnn- 2 ) and that of DSF is 0.92 Wnn- 2 (standard 
deviationO.57 Wnn- 2 ). The mean total OLR in the second case (column D) differs by -0.39 
Wnn- 2 (standard deviation is 1.36 Wnn- 2 ) and DSF differs by 0.86 Wnn- 2 (standard 
deviation is 0.62 Wnn* 2 ) compared to the line by line calculations. These values arc 
marginally poorer than those obtained by usingour model (column B in Table 3). Calculation 
CtnTahle3 shows that use of a fiequency-dqiendrait value df die zmife angle itssdf is not very 
critical. Calculation D in Table 3 shows that integrating monochromatic transmittances ova* 
angles and using these to-obtain a rapid algorithm is almost as good as finding band 
t ransmittanc es at effective angles which match the band radiative flUxes. 

.vi In addition to the above cases, we have carried out a calculation in which transmittances 
at nadir are averaged to get ban d transmittances, and then the rapid algorithm coefficients arc 
ftw H When calculating OLR and DSF, these band transmittances are then raised to the factor 
1£6. As shown in column E of Table 3, dearly, thlr approach Is very poor and introduces 
large biases in OLR and DSF. The bias is negative (-14.24 Wnn- 2 ) in QLR with a standard 
deviation of 2.78 Wnn- 2 , and positive (8.03 Wnn- 2 ) in DSF with a standard deviation of 
3.40 W na^ 2 , suggesting that absorption is overestimated in both calculations. This is expected 
because i*® 00 is always greater than , where the bar represents frequency average, 
except for the special case in which t is independent of v in the band. This last calculation of 
raising the band transmittances to 1.66, although not commonly used in band models, is 
out to indicate what would be the magnitude of error if thediffusivity approximation is 
mistakenly applied directly to bond transmittances instead of optical paths (effective absorber 

arhount). From the results presented in Table 3, it appears that if the angular integration is 

carried out for the monochromatic transmittances, and tile band parameters arc found from the 
angle-integrated transmittances, the method of exact angular integration is not critical. But, in 
the context of our band model, vtiudi does aot dhectly use optied mass, tiie approach we have 
used is correct and accurate. The benefits of this approach perhaps would be more evident if 
used in calculations of outgoing flux from narrow spectral brads. 

3.2 Comparison with the ICRCCM Results 

In the previous section, we have shown that our modd-cateulations agree very well 
with our line by Recalculations for clear-sky rad cloudy conditions. We now compare our 
model calculations of OLR and DSF with those of ICRCCM results for clear-sky rad cloudy 
conditions. We have used the standard atmospheric profiles of temperature, water vapor, and 




■ 
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ozone provided by the ICRCCM committee for tropical^ Mid-latitude summer ami winter, 
subarctic summer and winter conditions as compiled by McClatchey et al. (1972). For clear- 
sky conditions, we have calculated OLR and DSF using our model for cases 25, 27, 29, 31, 
and 33 of the ICRCCM tests (see Table 3 of Ellingson et al., 1991). Calculations of three line 
by line models, selected from more than 30 different models included in the ICRCCM, have 
been chosen for the comparison with our model calculations. The results of the comparison are 
shown in Table 4a. As shown in Table 4a, our band trrtdd calculations of OLR and DSF are 
within die range of 0.3% to 1.4% compared to the line by line calculations, and compare with 
the other line by line calculations as closely as they do to themselves. The major differences 
result from different assumptions about line and continuum absorption characteristics. 

For the cloudy case, we have used the mi&latitude summer profile and have carried out 
four different calculations of OLRand DSF/ as provided by the ICRCCM test cases 
corresponding to different cloud opacities. It is important to note that in the Pathfinder Path A 
calculations of OLR and DSF, we use the effective cloud amount which is a product of cloud 
emissivity and cloud fraction as determined from die TOYS measurements (Susskind et al., 
1987). Therefore, we do not have to explicitly parameterize cloud emissivity in the model for 
our OLR and DSF calculations. However, to test our model, we have parameterized cloud 
emissivity as a function of liquid water path and effective cloud-drop size following Paltridge 
and Platt (1981), Stephens (1978), and Mehta and Smith (1997). The model calculations 
carried out using this parameterization are presented in Table 4b along with the ICRCCM 
results. For these cases, line by line calculations from ICRCCM were not available. 
Therefore, the ICRCCM results are presented in terms of mean, standard deviation, and range 
of OLR and DSF based on the different models (Table 1 1 in Ellingson et al., 1991) which were 
used in the ICRCCM comparison. We find that our OLR calculations are within 1% and DSF 
are within 0.6% of the ICRCCM means few these cases. 

4.0 Summary 

in this report, the radiative transfer model used to calculate OLR and DSF from the 
TO VS Pathfinder Path A retrieval products has been described. The model is computationally 
efficient and is being used to calculate global, twice-daily per satellite, multi-year OLR and 
DSF using the TOVS Pathfinder Path A retrievals, which include surface skin temperature, 
atmospheric temperature, moisture, and ozone profiles, cloud top pressure and temperature, 
and effective cloud fraction a£ (Susskind et al., 1997). 

We have derived band-dependent effective zenith angles which are used for the angular 
integration of radiances in the model. These angles are used to calculate monochromatic 
transmittances at the appropriate angle which are then averaged to find band transmittances. 
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This is different from the conventional diffusivity approximation in which band transmittances 
as a function of optical path are found at nadir and the optical path of absorbers are scaled by 
the secant of a frequency-independent angle of 52.96°. The form of the rapid transmittance 
algorithm used for modeling the band transmittance is accurate for modeling wide band 
transmittances (100 to 300 cm-1) as used m the present study, as the model calculations differ 
by less than 0.2% from our line by line calculations. Our results arc also within 1% range of 
the results provided by ICRGCM for clear-sky and cloudy conditions, with some of the 
differences from FCRCCM resulting from effects of different assumptions of line and 
continuum absorption characteristics. 

' As shown by Mehta md Sussldad (1999), the QLR andCOLR derived from TOVS 
Path A arc m excelient agreement wMi theERBE measurements. We intend to validate the 
DSF with available measurements and will report the results elsewhere. The TOVS Pathfinder 
Path A QLR data ase currently available for 14 years (1985 to 1998) and will be extended to 
cover the period between 1979to 1998 in oearfiiiur& ~ 

> More information about the Goddard TOVS Pathfinder Path A data set is available 

from: http://faster.gsfc.nasa.gov/srt.html 
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Figure 1 : Spectrum of : a) outgoing longwave flux at top of atmosphere, b) brightness 

r • 

( ! temperatures corresponding to the outgoing flux, c) downwelling longwave 

flux at surface, d) brightness temperatures corresponding to the 
} i downwelling flux. Unit for flux is mW*nr 2 «cnr 1 , and for brightness 

y temperature is K. 



Figure 2 : Variation of outgoing longwave radiances (top) and downwelling radiances 
(bottom) with zenith angles. The five zenith angles shown are selected out 
of eight angles used for angular integration of radiances by Gaussian 
qudrature. Units are mW»m* 2 *ster 1 *(cnr 1 )' 1 . 

Figure 3 : Spectrum of effective zenith angles for outgoing flux (top) and for down- 
welling flux (bottom) calculated according to Equation 2.8. Units are degree. 

Figure 4 : Comparison of line by line and model calculations of OLR (left) and DSF 
(right) for 22 different profiles. Triangles in the diagrams show flux values 
under clear-sky conditions and squares show values under cloudy conditions. 
Units are W«nr 2 . > 
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Table la : Effective zenith angles averaged over 22 profiles aid their standard 
deviations for each model band used in the OLR calculations. 

Table lb : Same as in Table la but for DSF. 

Table 2a : Comparison of total OLR and its spectral components, averaged over 22 

profiles ami their standard deviations, with those from line by line calculations 
undo' dear conditions. 

Table 2b : Same as in Table 2a but under cloudy conditions with unit cloud emissivity 
and cloud top placed at 500 mb level in each profile. 


Table 2c : Same as in Table 2a but for DSF. 

Table 2d : Same as in Table 2b but for DSF. 

Table 3 : Comparison of total OLR and DSF, averaged over 22 profiles and their 
standard deviations, calculated by using different methods of angular 
integration, with those from line by line calculations, in which radiances were 
integrated over eight zenith angles by Gaussian qudrature. 

Table 4a : Comparison of model calculations of total OLR and DSF with those from 
ICRCCM project results under clear conditions. 

Table 4b : Comparison of model calculations of total OLR and DSF with those from 
ICRCCM project results under cloudy conditions. 
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Table la 


Equivalent Zenith Angles for OLR 


Band Equivalent 

(cm'l) Zenith Angle (St.Dev) 

degree 



0e (4®e) 

see (9 e ) 

2-150 

50.36 (0.59) 

1.57 

150-250 

50.64 (0.65) 

1.58 

250 - 350 

51.62 (0.53) 

1.61 

350 - 500 

52.57 (0.48) 

1.65 

500 - 650 

52.53 (0.37) 

1.64 

650 - 800 

52.94 (0.09) 

1.66 

800 - 950 

55.86 (0.78) 

1.78 

950 - 1 100 

53.08 (0.29) 

1.67 

1100-1250 

55.21 (0.26) 

1.75 

1250- 1500 

51.64 (0.33) 

1.61 

1500 - 1800 

50.36 (0.33) 

1.57 

1800-2100 

52.51 (0.41) 

1.64 

2100 - 2400 

52.74 (0.12) 

1.66 

2400 - 2750 

55.08 (0.14) 

1.69 


Table lb 



Equivalent Zenith Angles for DSF 


Band 

(cm**) 

EgiSvalent 

Zenith Angle (St J>ev) 
degree 
0e (A0«) 

sec (6 e ) 

2-100 

44.99 (0.36) 

1.41 

100 - 250 

44.54 (0.28) 

1.40 

250 - 350 

44.20 (0.33) 

1.39 

350-500 

43.79 (1.57) 

1.38 

500 - 650 

46.04 (2.24) 

1.44 

650 - 800 

52.12 (0.15) 

1.63 

800 - 950 

55.71 (0.23) 

1.77 

950-1100 

5432 (0.19) 

1.74 

1100- 1250 

53.92 (0.56) 

1.70 

1250 - 1500 

47.45 (1.14) 

1.48 

1500- 1800 

42.88 (0.74) 

1.37 

1800 - 2100 

50.57 (0.58) 

1.57 

2100 - 2400 

52.38 (0.17) 

1.64 

2400-2750 

54.85 (0.07) 

1.74 


Comparison of Model Outgoing Longwave Radiation with Line-by-Line Calculations 

Clear-sky Condition 
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Comparison of Model Outgoing Longwave Radiation with Line-by-Line Calculations 
Cloudy Condition : 100 % cloud cover with top at 500 mb. cloud emlssivity = 1.0 




Comparison of Model Downwelling Surface Longwave Radiation with Line-by-Line Calculations 

Clear-sky Condition 
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Comparison of Model Downwelling Surface Longwave Radiation with Line-by-Line Calculations 
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A) Une-by-line calculations, B) Model calculations as shown in Table 2a and 2c, C) Monochromatic transmittances are found at 52.96° and then 
the model coefficients are found for each band, D) Angular integration of monochromatic transmittances over all angles is carried out by Gaussian 
quadrature and then the model coefficients are found, E) Band transmittance coefficients are found at nadir and then the transmittances are raised to 
the factor 1.66. 



Table 4 li ■ . 

Comparison of OLR and DSF with ICRCCM Clear-sky Results 


Case 

OLR 

(W*m-2) 

OLR 

(W-m-2) 

DSF 

(Wm- 2 ) 

DSF 

(W.m- 2 ) 


Present 

Study 

From 

ICRCCM 

Present 

Study 

From 

ICRCCM 

Tropical 

296.69 

294.48 a 

298.28 b 

394.46 

389.89* 

389.46 b 



296.83° 


389.96° 

Midlatitude Summer 

289.89 

’ 285.78* 

288.99 b 

344.84 

342.31* 

340.20 b 



288.05° 


343.23° 

j Midlatitude 
Winter 

233.05 

234.26* 

236.56^ 

217.78 

215.42* 

210.54 b 



235.23° 


217.27° 

Subarctic 

Summer 

266.52 

266.29* 

269.96 b 

294.45 

292.79* 
287.3 l b 



269.44° 


293.42° 

Subarctic 

winter 

200.16 

200.87* 

202.95 b 

165.24 

164.45* 

162.45 b 



201.27° 


168.23° 

a Line-by-line model - Schwarzkopf and Fels 

b Line-by-line model- Scott and Gnedin 

c Line-by-line model - AiidngA. 




Table 4b 


Comparison of OLR and DSF with ICRCCM Cloudy Results 


Case 

OLR 

(W.m- 2 ) 

OLR 

(W*m- 2 ) 

DSF 

(W.m- 2 ) 

DSF 

(W.m- 2 ) 


Present 

Study 

From 

ICRCCM 

Present 

Study 

From 

ICRCCM 

Cloud top at 13 km 
liquid water 10 g*m 
effective dropsize 5.25 Jim 

163.76 

165.06 a 
18.75 b 
66.28 c 

359.05 

360.43 a 

6.74 b 

23.10° 

Cloud top at 2 km 

2 

liquid water 10 g«m 
effective dropsize 5.25 |im 

275.19 

276.05 a 

8.38 b 

33.68 c 

400.23 

398 .97 a 
8.02 b 
28.08° 

Cloud top at 13 km 

2 

liquid water 200 g*m 
effective dropsize 31 |im 

128.35 

129.82 a 

5.26 b 

21.98 c 

361.16 

361. 10 a 
7.22 b 
24.22° 

Cloud top at 2 km 

2 

liquid water 200 g*m 
effective dropsize 5.25 Jim 

271.02 

273.47a 

8.68 b 

39.67° 

410.18 

412.69 a 

2.54 b 

7.91° 


a 

b 


c 


Mean value from die models included in the ICRCCM tests 
Standard Deviation from the mean 
Range of values found in the models 
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